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The bb production cross section in 920 GeV proton - nucleus collisions was obtained with the 
HERA-B detector, using a sample of ~ 1.35 million of di-lepton triggered events taken during 
the 2000 physics run. The bb events are identified via the inclusive decay b — > J/ip, where 
the J/t{j mesons decay to fj>~ or e + e~ . The b events are separated from directly produced 
J/i/>'s by requiring that the J/ip decay vertex does not coincide with the primary vertex. The 
total measured cross section per nucleon,combining the results obtained in the two channels, 
is a(bb) = 32^ 1 2(stat)^(sys) nb/nucleonM The result is compared with the most recent QCD 
prediction and other experimental results. A much more precise measurement will be obtained 
from the 2002/2003 data. 



1 Introduction 

The heavy quark production cross section provides a good test of perturbative QCD theory. 
Theoretical calculations have been obtained up to the next to leading logarithmic order (NLLO), 
including the re-summation of soft-gluon effects^ The uncertainties of these results are still 
quite large. Several experimental measurements have been performed to test production and 
hadronization of the heavy quarks. Although most of these results are in qualitative agreement 
with theoretical predictions, there are some quantitative discrepancies that may be solved with a 
better understanding of non-perturbative processes in heavy quark production and hadronization 
mechanisms or introducing new physics effects^. 

Experimentally, the hadroproduction of the beauty quarks has been studied by several collab- 
orations at CERN and Fermilab. For proton - nucleus collisions there are only two measurements 
of the bb production cross section, done by the E789 and E771 collaborations. E789 has obtained 
a b i = 5.7 ± 1.5(stat) ± 1.3(sys) nb/nucleon, for 800 GeV p- Au collisions, using the decay chain 
b — * J/ipX — > fj, + fj,~X to identify the bb pairs^J. At the same proton beam energy, E771 has 
measured a h i = 43 ± 7 nb/ nucleon by detecting the two muons coming from the semileptonic 
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Figure 1: Layout of the HERA-B detector. 



decays of the b and b quarks and using Si as a target^. Although these results are compatible 
within the large experimental errors, they differ by around a factor seven. HERA-B can provide 
another measurement for 920 GeV proton - nucleus collisions, reducing the experimental errors 
in this energy region. 



2 The HERA-B experiment 

HERA-B is a forward spectrometer with a wide angular coverage, designed to measure CP 
Violation in the B mesons system^ Using a fixed target in the halo of the 920 GeV proton 
beam of the HERA storage ring (DESY, Hamburg), the spectrometer was designed to cope with 
very high interaction rates. It is formed by a vertex detector, a tracking system and particle 
identification detectors (see figure P). The target system consist of two stations made of four 
wires each, separated by 4 cm along the beam direction. The wires are made of different materials 
and can be individually moved to the halo of the proton beam untill the desired interaction rate 
is reached. 

The Vertex Detector System (VDS) was designed to accurately separate the decay point of 
B mesons from the main interaction point. It is formed by eight stations of four double sided 
silicon micro-strip detector planes (four stereo views), with a pitch of 50 /im. The first seven 
stations are in a Roman Pot system, inside the beam pipe, and can be displaced closer to the 
beam during the running period. The vertex detector reached design performance during the 
2000 run, with a J ftp vertex resolution of 60 /um in the direction transverse to the beam and 
500 //m in the longitudinal direction. 

Behind the VDS, the tracking system is placed. It is formed by 13 tracking stations that use 
two different technologies for the inner region, where the particle flux is larger, and for the outer 
region. The Inner Tracker (ITR) uses a micro-strip gaseous detector with GEM (Gas Electron 
Multiplier). It was in a commissioning phase during the 2000 data taking period, so it was not 
used in the trigger and consequently was also not used in the analysis presented here. The Outer 
Tracker detector (OTR) is made of honeycomb drift chambers, with a wire pitch of 5 mm in the 
innermost part and a 10 mm pitch in the outermost part. During the 2000 run, it reached a hit 
efficiency of 90% and a hit resolution of 350 |jm. 

The particle identification is provided by three different subdetectors: the ring imaging 



Cerenkov detector (RICH), the electromagnetic calorimeter (ECAL) and the muon chambers 
(MUON). The RICH can provide a 4a separation for e — n between 3.4-15 GeV, tt — K between 
10-60 GeV and K — p between 23-85 GeV. The ECAL, placed after the tracking system, provides 
electron and photon identification. It has a shashlik structure with Pb and W as absorbers, 
sandwiched between scintillator layers. The measured energy resolution during the 2000 run in 
the inner part was a{E)/E = 22.5%/E © 1.7% in the range from 12 to 60 GeV, close to the 
design values. The spatial resolution was 0.2 cm. 

Finally, the muons are identified by the muon system, made of four stations of tracking 
detectors placed in between iron absorbers. Different technologies are used, depending on the 
distance to the beam center. The inner part is covered with gas pixel chambers, while the outer 
part uses conventional tube chambers. In addition to this, there were pad chambers in the outer 
part that were used at the trigger level. 

The HERA-B trigger system was designed to identify the J/ip — > l + l~ decay coming from 
the decay of the B mesons, in both the e + e~ and the channels. The trigger chain is started 

by a signal coming from the muon system or the electromagnetic calorimeter (pre-trigger signal) , 
that may be a high pr track (two pad coincidence) in the muon system or a high transverse 
energy cluster in the ECAL. The First Level Trigger (FLT) uses the signal from the pre-trigger 
to start a fast tracking through the tracking chambers, based on regions of interest. The Second 
Level Trigger (SLT) confirms the FLT tracks by applying a refined tracking algorithm, based on 
the Kalman filter, and continuing the tracking to the vertex detector. Once the event is accepted 
it is fully reconstructed online and classified according to different interesting analyses. 

During the 2000 run the FLT was still in a commissioning phase, so it was not actually 
filtering and the SLT had a larger reduction factor than originally designed. The trigger system 
was able to reduce the rate from 5 MHz to 20 Hz. In this condition, ~ 0.9 M di-muon and 
~0.45 M di-electron triggers were collected during a short physics run, using two different target 
materials: C (77% of the data) and Ti (23%). 



3 Method 



The procedure used to measure the bb cross section consists of taking advantage of the long 
lifetime of the B mesons (~ 8mm in the HERA-B conditions). The idea is to identify the 
reaction pA — > bbX — > J/ipX, where the J/ip decays to /U + /x~ or to e + e~, by requiring that 
the J/ip decay vertex should be clearly separated from the primary interaction point. The good 
vertex resolution is, therefore, essential for this analysis. 

The bb cross section is measured relative to the prompt J/ip production cross section in 
order to reduce systematic errors due to luminosity and trigger efficiency. The acceptance of the 
HERA-B detector in xf is mainly in the negative range, from -0.25 to 0.15, in contrast to the 
other experiments which were only sensitive to the forward region. The bb cross section in the 
HERA-B acceptance can be writen as: 

A < N B 1 

Aaj /ip N P e R e^Br(bb - J/iPX) 1 ' 

where cx^ and Oj/^, are the b — > J/ip and the prompt J/ip production cross sections in the 
HERA-B acceptance, for a target with atomic mass A; Nb and Np are the number of detached 
and prompt J/ip decays, respectively; ep is the ratio between the reconstruction efficiency for the 
detached and the prompt J/ip decays; e^ 2 is the efficiency of the detached selection requirements 
and Br(bb —* J/ipX) is the branching ratio of the bb pairs decaying into a J/tp plus something 
else, that was assumed to be the same as measured for the Z decays: 2 • (1.16 ± 0.10)%^. 

The prompt J/ip cross section was measured by E789 for p-Au collisions'*^ and by E771 
for p-Si collisions . The value used in this analysis was obtained by re-scaling the Fermilab 
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Figure 2: /-i + fi~ invariant mass spectrum after the J/ip selection requirements. The solid line shows the fit obtained 
with the sum of a Gaussian for the signal and an exponential for background. The dashed line represents the 
background obtained with like-sign combinations that differs from the background under the J/ip peak due to 
the different trigger acceptance for the like/unlike-sign combinations and due to the physics contributions to the 

unlike-sign spectrum (Drell-Yan, open charm). 

measurements to the HERA-B proto n b eam energy^^, correcting for the most recent atomic 
mass dependence (a = 0.955 ± 0.005)QH anc l 

averaging the two results. The final prompt J/ip 
cross section is: (ij/^ = aj^/A a = (357 ± 8(stat) ± 27(sys)) nb/nucleon. Since for open charm 
and beauty production no nuclear suppression is expected and in the case of open charm it was 
not observed!*^, a is assumed to be 1 for the bb production cross section. 

The efficiencies of the reconstruction and detached vertex selection are calculated with the 
Monte Carlo simulation. The HERA-B MC uses Pythia^^to generate the production and the 
subsequent hadronization of the heavy quark pair in the hard interaction of the proton with 
the target nucleus (pA — > QQX). The remaining interactions inside the nucleus are generated 
with the package Fritiof^l To properly simulate the J/ip kinematics, the distributions of the 
differential J/ip cross sections as a function of if and pr, measured in p-Au collisions 
were used to weight the generated events. For the bb production a recent theoretical model, 
based on the calculations by Mangano et aZ. 021 and 

using the most recent NNLL MRST parton 
distribution functions ^!, was used to weight the generated events. Then, the generated particles 
are passed to GEANT which simulates the HERA-B detector, including mapping of the dead 
channels and realistic hit efficiencies. A full trigger simulation is applied and finally, the MC 
data is reconstructed with the same procedure used for the real data. 

The theoretical parameters in the MC simulation were varied within the range predicted by 
theoretical uncertainties. The variations induced in the bb cross section were taken a systematic 
error and included in the final measurement^. 

4 Prompt J /ip selection 

In order to identify prompt J/ip decays, the trigger tracks are used offline to seed the track 
reconstruction algorithm in the OTR. The offline-reconstructed tracks, including the segment 
inside the VDS detector, are then matched to the trigger tracks using a x 2 criterion. The 
selection requirements for the J/ip decaying into the electron channel are different from those 
applied to the muon channel due to the different triggering conditions and background levels. 

For the muon channel, additional particle identification requirements are applied to the 
offline reconstructed trigger tracks, taking into account the information provided by the muon 
chambers and the RICH. A good quality secondary vertex is also required. The resulting J/ip 
mass spectrum is shown in figure [5J 2880 ± 80 prompt J/ip — > decays were reconstructed. 

For the electron channel, the background is much larger so stronger particle identification 
criteria are used. The first one requires that the energy measured in the ECAL should be 
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Figure 3: e + e~ invariant mass spectrum for different requirements on the bremsstrahlung photons: no special 
requirements (left), at least one bremsstrahlung photon emitted (middle) and two bremsstrahlung photons emitted 
(right). In all the cases, the solid line represents the fit to a Gaussian signal plus a polynomial background. The 
dashed lines represent the background due to a wrong photon association. 



equal to the momentum measured by the tracking chambers (\E/p — 1| < 3a, with a ~ 9%), 
as the electrons are expected to leave all their energy in the ECAL. The second requirement 
identifies the electrons by the emission of a bremsstrahlung photon. When the bremsstrahlung 
photon is emitted in the VDS, it keeps the original direction of the electron track untill it is 
stopped at the electromagnetic calorimeter, providing a clear signature that can be used to 
isolate the J/tp — ► e + e~ signal. The power of this requirement is illustrated in figure 01 where 
the e + e~ invariant mass was plotted for three cases: when no bremsstrahlung photon is required 
(left), when one of the electron tracks is required to have emitted a bremsstrahlung in the VDS 
(middle) and when both electron tracks are required to have emitted a bremsstrahlung photon 
in the VDS (right). The efficiency of the bremsstrahlung reconstruction was measured using 
the plots from figure El and found in agreement with the MC simulation. The obtained value is 
(-brem = 0.34 zt 0.02(stat) ± 0.02(sys). tbrem can be used to infer the number of prompt J/ij) — ► 
e + e~ decays when no bremsstrahlung requirement is applied: Np = 5710±380(stat) ±280(sys). 



5 Detached J/ip selection 

The b — > J/ip events are selected using the long decay distance of the beauty hadrons, that imply 
also large impact parameters of the daughter tracks with respect to the primary interaction point. 
As the beauty hadrons are mainly produced along the beam direction, the longitudinal decay 
distance is used. The optimization of the selection requirements is done by maximizing the 
significance of the signal, S/y/B, where S is the number of signal entries from the MC and B is 
the background from the real data. 

For the muon analysis, the longitudinal decay distance between the secondary vertex and 
the target wire is required to be larger than 7.5 cr; the impact parameter of the muon tracks with 
respect to the wire, I w , should be larger than 45 //m and the impact parameters of the muons 
with respect to the primary vertex, I pr im, should be larger than 160 fim. Both requirements on 
Iprim and I w are needed, the first one because a 2-dimensional impact parameter gives better 
separation power and the second one to reject di-muon candidates that may have been produced 
at other primary interaction vertices on the same wire. 

For the electron analysis, the longitudinal decay distance calculated with respect to the 
target wire is required to be larger than 0.5 cm. The impact parameters of the electron tracks 
with respect to the wire should be larger than 200 fim or the minimum distance between the 
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Figure 4: (upper) and e + e~ (lower) invariant mass spectrum after the detached selection requirements. 

The upstream (unphysical) region is shown on the left and the downstream spectrum on the right. The solid line 

represents the unbinned likelihood fit (see text). 



electron track and any other track at the position of the wire should be larger than 250 /im. 
The isolation requirement plays the same role as the cut on the maximum impact parameter 
with respect to the primary vertex, while keeping the efficiency higher. The efficiencies of the 
detached selection conditions and the relative reconstruction efficiencies, as obtained with the 
MC, are e^e^ 2 = 0.41 ± 0.01 for the muon channel and €r€^ z = 0.44 ± 0.02 for the electron 
channel. 

The invariant mass spectrum for the ^ + [T~ and e + e~ pairs are shown in figure 0] The com- 
binatorial background, simulated by the events found in the region upstream from the primary 
vertex (unphysical region), is also shown there. The number of signal events was obtained with 
an unbinned maximum likelihood fit. In the case of the muon channel, a Gaussian function for 
the signal plus an exponential background were used. The parameters of the Gaussian were 
taken from the prompt J/ip spectrum. The resulting number of b — > J/ip — ► f^ + fJ-~ decays is 
1.9^'g. For the electron channel, the signal shape was taken from the MC b — > J /if), while the 
background was a combination of the shape from the upstream spectrum and the MC double 
semileptonic b decays. The results of the fit yields 8.6^32 b — > J/ifj decays. Using the efficiencies 
shown before and after a weighted average of the two target materials, the resulting bb cross 
sections in the HERA-B acceptance region are Ao" fe g = 16[ti2 nb/nucleon for the muon channel 
and Acj fe 5 = 38^^ nb/nucleon for the electron channel. The two results are compatible within 
statistical uncertainties. 

The results obtained in both channels were combined performing a common maximum like- 
lihood fit with four parameters (Acr^, A*" 1 "/" - background slope, e + e~ and n + fi~ yields). The 
result from the fit is Aer 6 g = 30^5 (stat) nb/nucleon. The main systematic uncertainties of this 
result are due to the prompt J ftp cross section measurement (9%), the trigger and detector sim- 
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Figure 5: Comparison of the bb production cross section measured by the HERA-B collaboration with the theo- 
retical predictions and the other experimental results. 
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Figure 6: Preliminary detached peaks obtained for the J/tjj — > e e (left) and J/tjj 
channels after processing par of the data from the 2002/2003 run 
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ulation (5%), the prompt MC production models (2.5%), the bb MC production models (5%), 
the prompt J /ip — * e + e~ counting and the difference in the C and Ti efficiencies (1.7%). The 
uncertainty due to cut values and background shapes in the fit are 13% for the muon channel 
and ^24% for the electron channel. 

Extrapolating to the full xp region, the total bb production cross section for 920 GeV proton - 
nucleus collisions is: 

a bl = 32±i!(stat)lf(sys) nb/nucleon (2) 

This result is compared with the theoretical predictions at the next to leading logarithmic order 
and with the other experimental results in the figure The HERA-B measurement is in good 
agreement with the theoretical prediction and is compatible with the previous experimental 
results. 



6 Prospects for the new data 

HERA-B has continued taking data during the year 2002 and beginning of 2003. A total sample 
of ~170 M di-lepton triggers were acquired, corresponding to an expected yield of ~ 300 K prompt 
J /if) decays. With this data a new measurement of the bb production cross section is expected to 
come soon. Preliminary detached J/tj) peaks were obtained already for both J/ip decay channels, 
after processing part of the data (see figure EJ) • 



7 Conclusions 



A sample of ~ 1.35 M di-lepton triggers taken during a short physics run in the year 2000 
were analysed. By exploiting the long lifetime of the beauty hadrons, a total yield of 8. 6jl|;| 
6 — > J/ij) — ► e + e~ and 1.9li'§ ^ ~~ * ■W' ^ / ti+ / u_ decays were obtained. The 66 production cross 
section was measured relative to the prompt J/tp cross section. The combined result for the 
two channels, after extrapolating to all the xf range is a b i = 32±;[2( s tat)i7(sys) nb/nucleon, in 
good agreement with the next to leading logarithmic order QCD predictions. A more accurate 
measurement is expected to be obtained soon, after analysing the data acquired during the 
2002/2003 run. 
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